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Abstract 

We investigate solutions to the equation dt£ — V/S.E = XS'^S, where S{x,t) 
is a Gaussian stochastic held with covariance C{x — x',t,t'), and x S W^. 
It is shown that the coupling XcN{t) at which the A^-th moment {£'^{x,t)) 
diverges at time t, is always less or equal for T> > than for D = 0. Equality 
holds under some reasonable assumptions on C and, in this case, XcNit) = 
NXc{t) where Xc{t) is the value of A at which (exp[A Jq S'^{0, s)ds]) diverges. 
The V = case is solved for a class of S. The dependence of XcN{t) on d 
is analyzed. Similar behavior is conjectured when diffusion is replaced by 
diffraction, D — > iD, the case of interest for backscattering instabilities in 
laser-plasma interaction. 
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I. INTRODUCTION 



We investigate the development of a linear amplification in a system driven by the square 
of a Gaussian noise. This problem arose and continues to be of interest in modeling the 
backscattering of an incoherent high intensity laser light by a plasma. There is a large 
litterature on this topic, and we refer the interested reader to Ref. |jT[| for background. 
Our starting point here is the work by Rose and Dubois |0 who investigated the following 
equation for the complex amplitude S{x, z) of the scattered electric field 

d,£{x,z) -iV/\£{x,z) = X\S{x,z)\^S{x,z), 
ze[0,L], a; G A C M^ and S{x, 0) = So{x). 

In Eq. (|I]), z and x correspond to the axial and transverse directions in a plasma of length 
L and cross-sectional domain A. The input at z = 0, Sq{x), is a given function of x and A 
will be generally taken to be a torus (e.g. in numerical solutions of Eq. ([l|) using spectral 
methods). The coupling constant A > is proportional to the average laser intensity and 
P is a constant parameter introduced for convenience. The complex amplitude of the laser 
electric field S{x, z) is a homogeneous Gaussian stochastic field defined by 

{S{x,z)) = {S{x,z)S{x',z')) = 0, 
{S{x, z)S{x', z'Y) = C{x -x',z- z'), 

where the correlation function C{x,z) is the solution to 

d,C{x,z) + ^AC{x,z)=0, 

(2) 

z e[0,L], X e A, and C{x, 0) = C{x), 

with C{x) a given function of x normalized so that C(0) = {\S{x,z)\'^) = 1. 

Using heuristic arguments and numerical simulations, Rose and DuBois found that the 
expected value of the energy density of the scattered field {\S{x,L)\'^) diverged for every 
L > as A increased to some critical value Xc{L). The average {\S\'^) is over the realizations 
of the Gaussian field S. This divergence indicates a breakdown in the assumptions made in 
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deriving Eq. (|T|), which neglects both nonhnear saturation and transient time evolution 0,0. 
Physically, it can be interpreted as indicating a change in the nature of the amplification 
caused by the plasma. 

To see the origin of this divergence in its simplest form, consider the case where V is set 
equal to zero in Eq. (p, and neglect all dependence of S' on x and z. We are then led to the 
equation 

which yields 

S{z) = £:(0)e^^'^ z>^. 

Here 5*^ = 5*^ + 5*1 and 5*1, 5*2 are two independent real Gaussian random variables with zero 
mean and unit variance. It is easily seen that the probability distribution of £{z), setting 
£^(0) = 1, has the density 

W{8,z) = (2A2)-^^~[^+(2^")"'l for 8>l, z> 0. (4) 

If we now take moments of £ at some value L of we find that {S^{L)) will diverge 
whenever 2NXL > 1. At the critical coupling Xcn{L) = {2NL)~^, there is a qualitative 
transition of the amplification of {S^{L)) from a regime where it is dominated by the bulk 
of the order-one-fiuctuations of 5* to a regime where it is dominated by the large fiuctuations 
of S. This toy model can be thought of as an idealization in which the size of the plasma 
is very small compared to the correlation length of the laser field. This is certainly not a 
reasonable physical approximation and we shall later consider situations in which S in Eq. 
(0) is 2;-dependent with a covariance C{z, z'). The equation is then still solvable more or less 
explicitly, depending on the form of C, at least as far as the dependence of the divergence 
of the moments of £^ on A and L is concerned. The main difference from Eq. is that for 
small enough values of A, the first few moments need not diverge for any L. 

In this paper, we extend these results to the x-dependent case where iT> in Eq. (|T]) is 
replaced by V, i.e. we consider a diffusive process in x rather than a diffractive one. Some- 
what surprisingly the diffusion does not suppress the onset of divergences in moments of the 
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field. This suggests a similar behavior for the diffractive case - in accord with the numerical 
results of - but we are unable to prove this at the present time. 

Before going on to the formulation and presentation of results for the diffusive case, we 
make some remarks about the relation between expectations over different realizations of 
the Gaussian driving term IS'p and the outcome of a given experiment. Accepting the 
idealizations inherent in assuming Gaussian statistics and neglect of nonlinear terms, the 
physically relevant question relating to the solution of the stochastic PDE (||) appears to be 
the following: What is the probability that for given A and L there will be small regions in 
A through which a significant fraction of the total incoming power is backscattered, (here 
"total" means through the whole domain A). Put more physically, imagine A to be divided 
up into M ^ 1 cells of equal area |A|/M and let R ^ 1/M be a specified number. We want 
to compute the probability P that in at least one of the cells the integral of \£\^ over that 
cell exceeds R\A.\. In the case where V is set equal to zero, this can be answered by taking 
for the cell size the transverse correlation length of \S\'^ and assuming the field inside each 
cell to be transversally constant and evolving along z under Eq. (j^) with a 2;-dependent S. 
One finds that P greatly increases as A passes its critical value for the divergence of the 
second moment, from P ^ 1 for A < \c2{.L) to P ~ 1 for A > \c2{L). We expect that this 
probability will behave similarly in real systems. 

The outline of the rest of this paper is as follows. In Sec. || we introduce our diffusion- 



amplification model. In Sec. |T| we prove that the value of the critical coupling obtained 
without the diffusion term cannot be less than the one obtained with the diffusion term. In 



Sec. |V| we prove that for a large class of Gaussian fields S the values of the critical coupling 
obtained with or without the diffusion term are the same. Section is devoted to the 
explicit solution of the diffusion-free problem in the particular case where the on-axis field 
S'(0, z) is a linear functional of a Gauss-Markov process. Finally, in Sec. |VT| we study the 
dependence of the critical coupling on the space dimensionality in the case of a factorable 
correlation function C. 
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II. MODEL AND DEFINITIONS 



As explained in the introduction, we consider a modified version of tlie linear convective 
amplifier model obtained by replacing iV by V on the left-hand side of Eq. (|l]). Taking 
V = 1/2 without loss of generality, one is thus led to the problem 

dtS{x,t) - ^AS{x,t) = XS{x,tfS{x,t), 

(5) 

t G [0, T], X e M*^, and £{x, 0) = £o{x), 

where, following the usual notation used in diffusion problems, the time variable t (resp. T) 
plays the role of the axial space variable z (resp. L). In Eq. (^), we restrict ourselves to the 
cases where S{x, t) is a real homogeneous Gaussian field defined by 

(^(x,t))=0, 

{S{x,t)S{x\t')) = C{x-x',t,t'), 

with the normalization C(0,0,0) = (^(a:, 0)^) = 1, and we will take Sq{x) = 1 as an initial 
condition. Note that S{x, t) is not assumed to be stationary in t, and that the rest of our 
analysis is essentially unaffected if we replace M'' by a (i-dimensional torus. 

The critical coupling \cn{T) and its diffusion- free counterpart AcAr(T) are defined by 

Xcn{T) = inf{A > : (^(0, Tf) = +oo}, (6a) 
A,^(r) = inf{A > : (e^^/o^^(0'*)''^*) = +oo}, (6b) 

where (.) denotes the average over the realizations of S. For a given T > 0, Eqs. (^ gives 
the value of A at which {S{x,T)^) blows up with and without diffusion respectively. 

Finally, in order not to make the calculations too cumbersome, we will use in the following 
the compact notation 
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t = (n,t), 

N T 
n=l 

S{t) = S{Xn{t),t), 

C(S, t) = {S{s)S{t)) = C{Xn{s) - Xm{t), S, t), 

Co{s,t) = C{0,s,t), 

i^,^lj) = j ip{t)^{t)dt, 

with s,t ^ [0)^]) n,m & N {1 < n,m < N), and where the are given continuous 

paths on M'^. The covariance operators Tq and tcg, respectively acting on (p{t) G L'^{dt) and 
ip(t) G L'^{dt), are defined by 

{fcip){s) = J C{s,tMt)dt, 

iic,ip)is) = r CiO,s,t)if{t)dt. 
Jo 



III. COMPARISON OF Xcn{T) AND Xcn{T) 

In this section we prove that AcAr(T) < AcAr(T). We begin with two technical lemmas 
that will be useful in the following. 

Lemma 1: Suppose the covariance function C{x,t,t') is continuous. Let 
... > be the eigenvalues of the covariance operator Tq- Here, the superscript x{t) denotes 
the continuous paths x„(t), 1 < n < A^. Then (exp A / S'(t)^dt) < +oo if and only if 
A < (2/ii*-*^)~^, and in this case one has 

log(a^J^")"»)^-li:iog(l-2A^f' ) < ^>-ffC(OMdt ^ 

To show d^), consider the Hilbert space of the L'^{dt) functions ^p{n,t) = ip(t) with 
the scalar product {ip^ip). Since C(s,t) is continuous in (s,t), and therefore bounded in 
compact sets, we have that / / C(s, t)^(is(it < +oo. By 0, Theorem VI.23, it follows 
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that the covariance operator is compact (and self-adjoint) in L'^{dt). Therefore there is 
an orthonormal basis {^j}j>i such that Tq^j = jXj^^^^pj. Consider now the sequence of 
random variables Xj = {S, (pj). As linear functionals of the Gaussian field S*, the Xj's form 
a Gaussian sequence with (X,,) = and (XiXj) = {ipi,Tc'Pj) = P-^^^^^ij. The equality in Eq. 
(0) is then obtained straightforwardly from / S'^{t)dt = Y^j':^Xj and the simple Gaussian 
identity (e^'^'^'^ = [\ — 2\jil^^^^ ^ , for 2A/i^'-*'' < 1. The inequality in Eq. follows from 
- log(l -x) < x/{l - x) and the fact that Ei fJ'f^^ = I C{t, t) dt = N C(0, t, t) dt. 

In the following subsection, Lp{t) = ip{n, t) will denote a set of N test functions normalized 
such that [if, if) = J2n=i lo V^(^' ^)'^ — ^■ 

Lemma 2: Assume that for every T > one has lim^^o sup^ jg^^-r] \C{x, s, t)—C{0, s,t) \ = 
0. Then, We > 0, 36 > such that 



< e 



for every G Bs^t, I < n < N, where Bs^t is the set of continuous paths x{-) such that 
\x{t) \ < 6 for every t E [0,T]. 

The proof of this lemma is straightforward: from the uniform convergence condition on 
C{x,s,t) it follows that We > 0, 35 > such that |C(s, t) — Co(s,t)| < e for every e 
^^T, 1 < ^ < A^- Thus, W > 0, 35 > such that 



< e' l^j \^{s)\dsj < e'NT, 



for every G B^^t, 1 < n < N . It remains to take e' = e/{NT), which proves Lemma 2. 

We can now state the main result of this section. Namely, that one of the diffusion effects 
on the divergence of the moments of S{x, T) is a lowering (or, more exactly, a non-increasing) 
of the critical coupling. The rigorous formulation of this result can be stated as the following 
proposition. 
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Proposition 1: For every T > 0, if limi,.^o sup^ jgjg.T] s,t) — C{0,s,t)\ = 0, then 
Ac7v(T) < JMT). 

In order to prove this proposition, one writes the moments of S in terms of the Feynman- 
Kac formula 



{£{0,Tf) = ((exp 



A J S{tfdt 



(8) 

x{t) 



where (■)x(t) denotes a A^-fold Wiener integral over Brownian paths 1 < n < A^, 

each arriving at x = 0. Let A > AcAr(T), i.e. fii > {2NX)^^, where fii is the largest eigenvalue 
of the covariance operator ic^- Let be the normalized eigenfunction associated with 

/ii, and (f){t) = (j){n,t) = N~^^'^(f)i{t) for every 1 < n < N. [N.B. : the factor N'""/"^ ensures 
the normalization (0,0) = 1]. By definition of one has 

/if^>(0,Tc0). (9) 

By Lemma 2, Ve > 0, 35 > such that 

(0, Tb0) > (0, Tbo0) - £ = iV/ii - £ (10) 

for every x„(-) G Bs^t, 1 < n < N. If one now takes e < Nfii — it follows from Eqs. (^ 
and (|10|) that Hi^^^ > 1/2A and so, by Lemma 1, 



(exp 



A J S{tfdt 



-oo 



for every x„(-) G Bs^t, I < n < N. Finally, since the set of the Brownian paths Xn{t) that are 
in Bs^T has a strictly positive Wiener measure, one finds from Eq. (H) that {S{0, T)^) = +oo, 
so A > Ac7v(T) which proves the proposition L 

Note that imposing the uniform convergence of C{x, s, t) to C(0, s, t) is not a very restric- 
tive condition. As far as we know, it seems to be fulfilled by any nonpathological stochastic 
field S of physical interest. 
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IV. EQUALITY OF Xcn{T) AND Acjv(r) FOR A CLASS OF S 



For a large class of Gaussian fields S it is possible to prove that diffusion has no effect 
on the onset of the divergence of {S{x,T)'^), i.e. Xcn{T) = Xcn{T). 

Proposition 2: Assume that for every T > one has linia^^osup^ j£[o,t] \C{x,s,t) — 
C{0,s,t)\ = 0, and that \C{x,s,t)\ < C{0,s,t) for every x e and s,t e [0,r]. Then 
Kn{T) — Xcn{T). 

The proof of this proposition is as follows: By the uniform convergence condition on 
C{x,s,t) and Proposition 1 one already has Xcn{T) < Xcn{T). It remains to show that 
Kn{T) < Xcn{T). Let III be the largest eigenvalue of the covariance operator tco- Let 0i(t) 
be a principal (normalized) eigenvector for the covariance operator Tq- One has 

/xf ) = (0i,Tb(/>i) < (|0i|, \fc<j>i\) < (|(/)i|,T'col0i|) < Nfxi, 

where the second inequality follows from the condition \C{x, s,t)\ < C(0, s, t). Suppose now 
A < Xcn{T), i.e. A < {2Nfii)-\ Then A < (2/ii'*V^ and by Lemma 1 



exp 



A J S{t fdt 



< cxp 



NXJoC{0,t,t) dt 



1 - 2A/it 



x{t) 



< exp 



iVA/o C(0,t,t) dr 
1 - 2NXni ^ 



Since this inequality is uniform over all Brownian paths, we finally have 

{S{0,T) )<exp(^ i-2NXi,i J ^ 

and therefore A < Xcn(T), which proves the proposition 2. 

This result shows that for Gaussian fields S fulfilling the not so restrictive conditions of 
Proposition 2, it is sufficient to solve the diffusion- free problem to determine the onset of the 
divergence of {£{x,T)^). It is therefore interesting to show how such fields can be actually 
obtained. To this end, the remaining of this section will be devoted to explicitely construct 
two typical examples of stochastic fields S which fulfill the conditions of Proposition 2. 
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(11) 



(12) 



A. an example of nonstationary S 

The first example is the diffusive counterpart of the Gaussian field defined by Eq. 
Let S{x,t) be the solution to 

dtS{x,t) - lAS{x,t) = 0, 
te[0,T], X e R'^, and S{x, 0) = S{x), 
where S{x) is a real homogeneous Gaussian field defined by 

(Six)) = 0, 

{S{x)S{x')) =C{x-x'), 
with C{x) a given function of x normalized such that C(0) = {S{x, 0)^) = 1. One has 

S{x, t) = j S{k)e'^''-"^''''d% (13) 

where S{k) is the Fourier transform of S{x), and from Eqs. (^) and ( |T3|) it follows that 
S{x, t) is a real homogeneous nonstationary Gaussian field with 

{S{x,t)) = 0, 

{S{x,t)S{x',t')) = /C(A;)e*'=(^-^'')-5^'(*+*')d'^A;, 
where C{k) is the Fourier transform of C{x). Since C{k) is real and positive 0, one has 

\Cix,s,t)\^\{Six,s)SiO,t))\ 
= I C(A;)e*^^-5fc'(^+*)rf'^A; 
< J C{k)e-'2''"('+'^d'^k = C{0,s,t), 

for every x G M*^ and s,t E [0, T], so S{x, t) fulfills the conditions of Proposition 2. 

B. an example of stationary S 



(14) 



The second example is provided by a modified version of Eq. (|ri|) obtained by adding a 
source term a la Langevin on its right-hand side. Namely, let S'(x, t) be the solution to 

dtS{x,t) - \/\S{x,t) = L{x,t), 
t e] - oo, T], X e R'^, and S{x, -oo) = 0, 



(15) 
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where the Langevin source term L{x,t) is a homogeneous Gaussian white noise defined by 

{L{x,t)) = 0, 

(16) 

{L{x, t)L{x', t')) = -6{t - t')A^C{x - x'), 
with C{x) a given ^ function of x normahzed such that C(0) = 1. The solution to Eq. ( p!5[ ) 
reads 

S{x,t) = I (fk y^"" f e-^^''(*-^)L(/t, s) dsl , (17) 

where L{k,t) is the Fourier transform of L{x,t). From Eqs. ([T6| ) and ([T7| ) it can be shown 
that S{x,t) is a real homogeneous stationary Gaussian field with 

{S{x,t)) = 0, 

{Six,t)Six',t')) = /C(fc)e''=("-"')-5'='M'lci'^fc, 
where C{k) is the Fourier transform of C{x). As previously, since C{k) is real and positive 
llSll, one has 



(18) 



\Cix,s,t)\^\{S{x,s)SiO,t))\ 

= J C{k)e'^''--^''^\'-'\<fk 

< J C{k)e-'^^"\'-'\d^k = C(0, s, t), 



for every a; G M'^ and s, t G [0,T], and so S{x,t) fulfills the conditions of Proposition 2. 
More generally, it can be checked that any real homogeneous Gaussian field S{x, t) defined 

by 

{S{x,t)) = 0, 

{S{x,t)S{x',t')) = S C{k,t,t'y^^^-'''U'^k, 
where C{k,t,t') is real and positive, fulfills the conditions of Proposition 2. 



V. EXPLICIT SOLUTION OF THE DIFFUSION-FREE PROBLEM FOR A 

CLASS OF 5 



In this section we show that an explicit computation of the diffusion-free amplification 
factor (exp(A^A Jq S{0,t)'^dt)) can be achieved if S{0,t) is a linear functional of a Gauss- 
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Markov process. Note that determining Xcn{T) amounts to finding the largest eigenvalue of 
the covariance operator , which in principle can always be achieved, at least numerically. 
As shown above, Xcn{T) > Xcn{T) with equality holding when Proposition 2 is applicable. 
Since Xcn{T) = N~^Xci{T) in the diffusion free case, we will take = 1 in the remaining 
of this section without loss of generality. 



A. Solution of the diffusion-free problem using the Feynman-Kac formula 

We consider the case where the Gaussian process S{0,t) can be written as 

S{0,t) = {c,Y{t)), (19) 

where (x, y) = x'^y = J2i ^iUi, c is a given n-dimensional vector, and Y{t) is a n-dimensional 
Gauss-Markov process defined as the solution to the linear stochastic differential equation 

dY{t) + AY{t)dt = GdB{t), 
F(0) Gaussian with (F(0)) = 0. 

Here, A and G are constant n x n matrices, and B{t) is a n-dimensional Brownian motion. 
From Eqs. (0) and it follows that one can write the diffusion-free amplification factor 
as a Feynman-Kac formula 

^^A/;5(o,t)^d*^ = ^e^/o"<^W'^^W)'^*^ =|t;(y,r)t/"i/, (21) 

where C denotes the symmetrical n x n matrix c (S> c, and v{y,t) is the solution to the 
parabolic equation 

f = {TtA + X{y, Cy))v + {Ay, V)v + \{G^V , G^V)v, ^^^^ 
^(?/'0) = (2^)"^'^exp [-\{y,U-'y)], 



I ■sn/2 ]^ 

^(?/'^)=(7^ ; exp 



with U = Cov[r(0), r(0)]. The solution to Eq. (^) has the form 

■^{y,K{tr'y) + xJ^TTCK{s)ds'^ , (23) 



27r 



\m\ 
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(24) 



where K{t) is a symmetrical n x n matrix which is the solution to 

= GGt - [AK{t) + K{t)A^ + 2XK{t)CK{t), 
K{0) = U. 
Thus, from Eqs. ( ^T]) and (l23| ) one has 

e^Io ^(0'*)' '^'^ = e^/o^ ^'^^^W (25) 

with -ft'(t) given by the Riccati equation (P^. 

The solution to Eq. ( P^ is known to explode in finite time for large enough A. For n = 1, 
in which case S{0,t) is itself Markovian, Eq. ( P^ is solved straightforwardly (see Sec. |V lj| ). 
For n > 2, the solution to Eq. ( p4|) can be obtained by the so-called Hamiltonian method: 
we define the 2n x 2n matrix 

^ At -2XC 
GG^ -A 



H 



/ 
V 



and solve the linear differential equation 

d 
It 

with the initial condition 



Qit) 


= H 


Q{t) 


_ P{t) 




_ Pit) _ 



(26) 



Q(o) 




/ 


_P(0)_ 




u 



The solution -ft'(t) to the Riccati equation (^^ is easily checked to be given by 

K{t) = P{t)Q{t)-\ (27) 

which explodes if and only if Q{t) becomes singular 0. Since Eq. (pB]) is a linear equation, 
it can in principle be solved by a symbolic computation program. 



B. Application to the n = 1 case 

As an example, let us consider the simplest case n = 1 with C{0,t,t') 
limit, the diffusion-free amplification factor reads 



e I* In this 



(28) 
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where Y{t) is the Ornstein-Uhlenbeck process 

dY{t) + Y{t)dt = V2dB{t), 



(29) 

(r(o)) = 0, {vm = 1, 



and K{t) is the solution to the Riccati equation 

1^ = 1- m + xm\ ^^^^ 

K{<d) = 1. 

Equation ( pOD can be easily solved by means of the substitution 2XK{t) = —d\ogu{t)/dt. 
Inserting the result into Eq. p8|), one obtains 



,xj:sio,t)-dt\ e^/^ ^^^^^^ ^^^^ 

'cosh(ar) + 2A) sinh(ar) 



.X s(o,tf dt\ e^/^ 

/ y'l + r/2 

/ y'cos(ar) + -2A)sin(aT) 
where a=\l- AX\^/^. It can be seen from Eq. (|53|) that, for A > 1/4, (exp(A S{0, tfdt)) 
diverges as T tends (from below) to the critical time Tc(A) given by 

where the determination of tan~^ is such that < tan^^ < tt. Inverting Eq. ([3^ ) and using 
\cn(T) = N^^\ci{T) gives the diffusion- free critical coupling Xcn(T) in the cases where 
CiO,t,t') =e-\'-''\. 

VI. DEPENDENCE OF THE CRITICAL COUPLING ON SPACE 

DIMENSIONALITY 

In this section we study the dependence of Xcn(T) on the space dimensionality D. We 
will restrict ourselves to the cases where the correlation function C can be written out as 

CDix,t,t') = Cd{x\\,t,t')CD-d{x±,t,t'), (35) 
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with Cd, Cd and Co-d continuous, symmetric, and positive definite, and where x\\ is the 
projection of x onto a rf-dimensional subspace {d < D) and x± = x — x\^. In the following, a 
correlation function of this type will be called a factorable correlation function. It is worth 
noting that such a correlation function can be very easily obtained, e.g. when the Gaussian 
field S is defined by either Eq. ( P^ or Eq. (^) in the cases where C{k) is factorable as 
C{k) = Cd{k\\)CD~d{k±). 

We prove that as A increases, the divergence of {£{x,T)^) obtained in the original D- 
dimensional problem cannot occur before the one obtained in the projected (i-dimensional 
problem whenever < (7^-^(0, t,t) < 1. Since many stochastic fields S of physical interest, 
e.g. in optics, do have a factorable correlation function, we expect this result to be useful 
for the comparison of two-dimensional numerical simulations with experiments and three- 
dimensional numerical simulations. Before expressing this result in a more rigorous way, we 
begin with two technical lemmas that will be needed in the following. 

Lemma 3: Consider a D-dimensional problem and let /ii*^*'* be the largest eigenvalue 
of the covariance operator T^^ and N given continuous paths x{t). Then Xcn{T,D) = 
[2 sup^(t)/if*^]"^ 

This lemma can be proven straightforwardly by successively considering the inequalities 
A > [2 sup^(t) f^i^^h~^ ^^'^ A < [2 sup^(t) A^i^*'*]"^) ^"^^ by following the same lines of reasoning 
as for the proofs of Propositions 1 and 2 respectively, where one replaces the N paths x{t) = 
corresponding to Xcn{T) = [2fil^^^~^]^^ by N paths maximizing fi^^^^ 0. 

Lemma 4: Let KQ{s,t), i^i(s,t), and i^2(s, t) be three symmetric kernels such that: 
(i) -^"0(8, t) = Ki(s, t)i^2(s, t); (ii) K2 is a positive definite continuous symmetric kernel; 
(iii) < K2{t,t) < 1 and the largest eigenvalue of Ki is positive, or K2(t,t) = 1 and no 
condition on the sign of the largest eigenvalue of Ki. Then fii{Ko) < fii{Ki), where fii{Ka) 
denotes the largest eigenvalue of Ka- 
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The proof of this lemma is as follows: since K2 is a positive definite continuous symmetric 
kernel, Mercer's theorem holds ||^ and this kernel admits the expansion 

ir2(s,t)=^a,/,(s)/,(t), (36) 

i 

where Oj > and /j(t) respectively denote the i^^ eigenvalue of the operator and the 
associated normalized eigenf unction. Let 0i(t) be a principal (normalized) eigenfunction of 
the operator Tk^ and ^i{Kq) the corresponding largest eigenvalue. From the condition (i) 
and Eq. (^), one has 

/ii(iro) = (01, TVo0i) = XI «i(/*0b TkiMi) = o-iMiiili, TK^rji), (37) 

i i 

where Mj and rii{t) are given by 

M, = (/A,/A), 

and 

77.(t) = M,-^/V.(t)0i(t), 

such that {rii,rii) = 1. By the definition of fii{Ki) and from K2{t,t) < 1, condition (iii), one 
has respectively 

m{Ki) > ivuTK,V^), (38) 



and 



ii(t)2rft 



= J K2{t,t)(f)i{tydt < J (j)i{tfdt = 1. (39) 

So, from Eqs. (pTf), (p^), (pOf ) and the condition (iii), it follows that iii{Kq) < fii{Ki), which 
proves Lemma 4. 

We can now proceed to rigorously express and prove the result stated at the beginning 
of this section. Let Xcn{T,D) be the critical coupling associated with a D-dimensional 
problem in which the correlation function of the Gaussian field S is given by C^i- One has 
the following proposition: 
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Proposition 3: for every T > 0, if CD{x,t,t') is a factorable correlation function such 

that < CD-d{0, t,t) <1 for <t <T, then Xcn{T, D) > Xcn{T, d). 

The proof of this proposition is straightforward. By the definition of a factorable corre- 
lation function one has Cd{^,'^) = Cd{s,t)CD-d{,^i^), where both Cd(s,t) and Co-di,^,'^) 
are continuous, symmetric, and positive definite. Since C£i_rf(t,t) = C/)_rf(0,t,t) and 
< CD_rf(0,t,t) < 1 by assumption, one can apply the lemma 4 with Kq = Cd, Ki = Cd, 
and K2 = Co-d- It follows immediately that /i^*-*'' < /ii where /ii denotes the largest 
eigenvalue of the operator T^^. Let Xmaxlt) be N paths maximizing jXl^^^ [^. One has 
suPa;(t) /Ui''*'' = /ii'"""^*^ < /ti'"'"'" from which it follows that sup^.^^) /ii''*'' < s\xpx^^(t)lH^^^^^ 
and, by Lemma 3, \cn(T,D) > \cN(T,d), which proves the proposition 3. 



VII. SUMMARY AND PERSPECTIVES 

In this paper, we have studied the effects of diffusion on the divergence of the moments 
of the solution to a linear amplifier driven by the square of a Gaussian field. We first proved 
that the divergence yielded by a diffusion-free calculation cannot occur at a smaller coupling 
constant than the one obtained from the full calculation (i.e. with diffusion). Then we have 
shown that, in the case where the absolute value of the (uniformly continuous) pump field 
correlation function is bounded from above by its one-point value, there is no diffusion effect 
on the onset of the divergence which is therefore given by a diffusion-free calculation. In 
this context, we have solved the diffusion- free problem explicitly when the pump field is 
a linear functional of a Gauss-Markov process. Finally, we have studied the dependence 
of the critical coupling on the space dimensionality in the case of a factorable correlation 
function. In particular, we have proved that the divergence obtained in a D-dimensional 
problem cannot occur at a smaller coupling constant than the one obtained in the projected 
(i-dimensional problem [d < D). 

As mentioned in the introduction, we would like to extend our results for the diffusion- 
amplification model (H) to the more difficult diffraction-amplification problem (|l]). According 
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to Eq. (H), the results obtained in this paper also apply, beside some minor technical modifi- 
cations, if the pump field is a complex Gaussian field as in Eq. (|1|). The remaining difficulty 
in extending our results to Eq. ([^) lies in controlling the complex Feynman path-integral, 
compared to that of the Feynman-Kac formula for the diffusive case. Expressing £{x, t) as a 
Feynman path-integral and averaging over the realizations of S*, one cannot a priori exclude 
the possibility that destructive interference effects between different path contributions make 
the sum of the divergent contributions finite. Thus one cannot deduce the divergence of the 
moments of S{x,L) from that of the amplification along paths arriving at the point {x,L). 
It is however not unreasonable to expect that Propositions 1, 2, and 3 also apply to the 
diffraction-amplification problem (|l|). Proving this conjecture is another matter and is the 
subject of a future work. Note that in the case of Proposition 2, the on-axis correlation 
function of the pump field must be real and positive, which is quite restrictive if the pump 
field is complex. From a practical point of view (e.g. in optics), it would therefore be very 
interesting to find out whether there exists an enlarged version of this proposition applying 
to complex on-axis correlation functions as well. 
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